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INTRODUCTION 


For many years, observations of the Earth radiation budget have been made from 
orbiting satellites. These observations can be analyzed by suitable techniques to 
estimate the radiant exitance field at the top of the atmosphere. A number of such 
data analysis techniques are available, each with its own simplifying assumptions and 
characteristics which affect the derived radiant exitance. In the past, there has 
been some consideration of how errors in measurements and errors associated with a 
given technique affect the derived radiant exitance. However, little emphasis has 
been placed on the differences between the techniques themselves and how these dif- 
ferences affect the derived radiant exitances. In the present study the effect on 
the radiant exitance estimates of analyzing the same data with three techniques is 
examined. 

The longwave component of the Earth radiation budget, measured by a wide field 
of view radiometer aboard an orbiting satellite, has been chosen as the data type for 
this study. Each measurement is an integral of the irradiance from all points within 
the field of view weighted by the directional response of the sensor. The shape 
factor technique has been the most popular way to solve this integral measurement 
equation. Each measurement is divided 1:^ a scaler to derive the radiant exitance at 
the top of the atmosphere. With another approach, the parameter estimation tech- 
nique, all measurements are processed together as a batch and the radiant exitance 
estimate is defined as a least sc^ares fit to the data. With the third approach, the 
deconvolution technique, use is made of the fact that spherical harmonics are the 
eigenfunctions of the integral measurement operator. These three techniques are 
studied by applying them to the same set of radiation data and comparing the result- 
ing radiant exitances on a global, zonal, and 10® regional scale. 

Details of the technique derivations and implementations are very important to 
understanding the results. For this reason the basic measurement equation which 
relates the measurements to the derived radiant exitance is formulated. Each tech- 
nique is derived from this equation, and the associated assumptions are set forth. 

In addition, each technicjue has its own unique implementation. For example, the 
results of the shape factor technique depend on the size of the surface area over 
which point estimates are averaged to produce regional estimates. The parameter 
estimation results are affected by the numerical integration scheme employed and by 
the model of the radiation directional function. The deconvolution results also are 
a function of this model and in addition vary with the degree of truncation of the 
spherical harmonic representation. Each of these aspects of the techniques is inves- 
tigated by means of numerical examples and discussed. The results of the deconvolu- 
tion technique are in terms of spherical harmonics, while the results of the other 
two techniques are regional radiant exitances. To compare these three techniques, 
their results must be transformed to the same coordinate system, and the transforma- 
tion also affects the results. Finally, we examine the advantages and disadvantages 
of each data analysis technique. 
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complex coefficients of spherical harmonics, W-m 
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S(a) 


cn' sn 
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area, km^ 

2 

area of globe, km 

matrix of influence coefficients (see eq. (13)) 

influence of radiant exitance of kth region on ith measurement 

integrals defined by equations (27) through (29) 

_2 

real coefficients of spherical harmonics, W-m 

shape factor (see eq. (7)) 

satellite altitude above R , km 

e 

total number of regions 

-2 -1 

Earth-emitted radiance, W-m -sr 
measurement operator 

—2 

measured radiation at satellite altitude, W-m 

-2 

radiant exitance , W-m 

degree of truncation; also, number of measurements 

normalizing coefficients for spherical harmonics 

associated Legendre polynomial of degree n and order m 

distance from surface element to satellite, km 

directional function, sr"^ 

radius of Earth-atmosphere system, km 

angular response of sensor 

arbitrary real function on the surface of a sphere 
real spherical harmonics of degree n and order m 
coitplex spherical harmonic of degree n and order m 

cone angle at satellite from satellite nadir to point on surface of Earth 
cone angle to horizon 
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P clock angle from north about the satellite nadir to point on surface 

of Earth 

smoothing parameter (see eq. (31)) 

Y Earth central angle 

Y* angular radius of a circular region 

6™ Kronecker delta function 

6 zenith angle 

0 colatitude 

nth eigenvalue of measurement operator 
p regional correlation coefficient 

p(Y) weighting function 

variance of random variable 

—2 

degree dispersion, W-m (see eq. (39)) 

4) azimuth 

^ longitude 

Q solid angle, sr 

Subscript: 
s satellite 

Abbreviations : 

FOV field of view 

LD limb darkening 

Lamb Lambertian 

A circumflex (^) over a symbol denotes an estimate. A tilde (~) over a symbol 
denotes the smoothed radiant exitance field. 


THEORY 

Formulation of Problem for Emitted Radiation 

The surface at the top of the Earth- atmosphere system is approximated by a 
sphere of radius R^. The Earth-emitted radiance L leaving any point on this 
spherical surface is modeled as a function of colatitude 0, longitude and zenith 


3 



angle 6 of the exiting ray and is not a function of azimuth (see fig. 1). The 

radiant exitance M at a point on this surface is given t)y 



Figure 1 Earth-satellite geometry. 


M(0,^) 




L(0,^,0) cos 0 sin 0 d0 d<() 


or 


r ^/2 

M(0,^) = 2n / L(0,^,0) cos 0 sin 0 d0 

*^ 0=0 


( 1 ) 


A radiation directional function for emitted radiation R(0/^/0)/ the limb darkening 
function, is defined so that 


L{0,$,9) = 1 M(0,®) R(0,®,0) 


( 2 ) 


In order for equations (1) and (2) to be conpatible, R must satisfy the normalizing 
condition. 


rV2 

2 / R(0,^,0) cos 0 sin 0 d0 = 1 (3) 

*^ 0=0 


For this study R is assumed to have no azimuthal dependence. If the surface is 
Lambertian and radiates with equal intensity in all directions, then R(0,^,0) = 1, 
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Now consider a radiation sensor at satellite altitude/ as shown in figure !• In 
modeling the satellite measurements, we do not consider measurement errors or time 
variation of the radiation field* Thus, the radiance from the top of the Earth- 
atmospheric system incident on the sensor at a colatitude Og, longitude and 

altitude h is modeled as 


m(0^,$ /h) = / L(0/^/O) s(a) dQ (4) 

s s ''FOV 


where 0 is a function of the locations of the satellite (0^,^^) and of the surface 
element at the top of the atmosphere (0,^)/ ^ is the solid angle at the satellite 

subtended by the surface element, a is the nadir angle at the satellite from the 
local vertical to the surface element, and the integration is carried out over the 
field of view (FOV) which depends on 0g, ^g, and h. The function S(a) is the 
angular response of the sensor to incoming radiation. This angular response function 
could be expanded to incorporate the geometry of other types of sensors as well as to 
incorporate a dependence on azimuth. However, in this study a perfectly black flat- 
plate sensor normal to the vertical is assumed for which S(a) = cos a. By use of 
equation (2), equation (4) may be written 


m(0 ,h) 

s s 


f M(0,$) R(0,$,e) cos a dQ 
■^FOV 


(5) 


which relates the known measurement m to the unknown radiant exitance M. 


Shape Factor Technique 

Many investigators (e.g. , Smith et al. 1977, Jacobowitz et al. 1979, Weaver and 
House 1979) have used the shape factor technique (also called inverse square law) to 
transform wide field of view radiation measurements at satellite altitude to esti- 
mates of the radiant exitance at the top of the atmosphere. The measurement m is 
divided by a geometric shape factor F to obtain the radiant exitance associated 
with the nadir point. The assumption is made that the radiant exitance at the top of 
the atmosphere is constant over the field of view, or that M(0,$) = M(0g,$g), so 
that equation (5) gives 


A(03,$3) 


m(0g,$g,h) 

F 


(6) 


where the circumflex ('') denotes an estimate and the geometric shape factor is given 
by 


^ i f R(©,®,9) 

Jpov 


COS oc dQ 


(7) 
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In most applications/ the shape factor is simplified hy the further assumption that 
the directional function is independent of position; that is, R(0,$,0) = R{0). If 
the variables of integration are changed from solid angle Q to the cone and clock 
angles a and P, then equation (7) becomes 


F 


2 

7C 




COS cz sin oc da dp 


or 


F 


2 


/““ R(e, 

a=o 


cos a sin a da 


( 8 ) 


where a^^ is the cone angle to the horizon, or the extent of the field of view. 
From figure 1 and the law of sines, the following relation is obtained: 


R 

sin a = ~ sin 6 (9) 

R + h 
e 


Differentiating equation (9) gives 

R 

cos a da = - — cos 0 d0 (10) 

R I 

e 


Substituting equations (9) and (10) into equation (8) yields 



or from equation (3) 
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Thus, the geometric shape factor is independent of the directional function for the 
assumptions considered and is a function only of the measurement altitude h. For a 
circular orbit, the shape factor is constant for all measurements. Whatever the 
case, the shape factor technique is characterized by reducing each measurement indi- 
vidually from satellite altitude to a radiant exitance at the top of the atmosphere 
hy dividing a geometric shape factor. These individual radiant exitances are used 
to define the global radiant exitance field. 


Parameter Estimation Technique 

With the parameter estimation technique, all measurements are processed as a 
batch, and a global radiant exit ant field at the top of the atmosphere is produced by 
means of a matrix inversion. Consider the surface of the globe divided into K 
regions, each with constant radiant exitance M, and directional function 
(k = 1, 2, . • • , K) . The measurement equation (5) can be e3q>ressed in terms of sur- 
face area instead of solid angle as 


m (0 ,h) = ^ f M( 0 , 3 i) r( 0 ,$, 0 ) cos a dA 

*^FOV r 


where r is the distance from the satellite to the surface element. Incorporating 
the assumed regional radiant exitance model gives 


m(0^,$^,h) = - 


k=1 -^FOV, 




cos a 


cos 


dA 


where the integration is over the surface area that is in both the field of view and 
the kth region. Thus, the ith measurement is modeled as 


K 


m 


i = ^ 

1 k=i ^ 


(i = 1, 2, N) (11) 


where is the influence coefficient of on the ith measurement and is given 

by 


. cos w . 

5 L„ “i 


ik 


dA 


( 12 ) 


^FOV., 

ik 


In matrix form, the set of equations defined by equation (11) become 


m = BM 


(13) 
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where m is a column vector of N measurements, ^ is a column vector of K 
regional radiant exitances, and B is the observational matrix with elements 

For a large measurement set (N > K) over the globe, equation (13) represents an 
overdetermined set of simultaneous equations which yield the least squares solution. 


^ T -1 T 
M = (B B) B m 


(14) 


Deconvolution Technique 

The deconvolution technique, described by Smith and Green (1981) and Bess et al. 
(1981), can be viewed as a parameter estimation technique with a spherical harmonic 
basis set rather than a piecewise constant basis set* The deconvolution solution to 
the measurement equation (5) can be found analytically, but it requires additional 
assumptions. It is convenient to express the relationship between the measured radi- 
ation unknown radiant exitance M(0,<^) as 


m(0 ) = ^[M(0,^] (15) 

s s 

where ^ denotes the linear integral measurement operator of equation (5). Smith 
and Green (1981) have shown* that the eigenfunctions of this linear operator are 
spherical harmonics, that is. 


se[Y”*(0,$)] = X Y®(0 ) 

n n n s s 


(16) 


where ^n^^s'^s^ ^ spherical harmonic of order 

the sxibsatellite point. The associated eigenvalue 


m and degree n 

A. is given by 
n 


evaluated at 


X 

n 


2 



P^(cos Y) R(0) cos a sin cc doc 
n 


(17) 


where P^(cos y) denotes the Legendre polynomial of degree n as a function of the 
Earth central angle y (see fig. 1). 

Because spherical harmonics are eigenfunctions of the measurement operator, let 
the radiant exitance at the surface be represented by a series of spherical harmonics 
truncated at degree N: 


M(0,$) = 2 y"*(0,S) 

n=0 m=-n 


( 18 ) 
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Also, let the measurements be represented as 


m(0,^) 


N 


E 

n=0 


E 


m=-n 


m 

a 

n 


n 


( 19 ) 


Substituting equations (18) and (19) into equation (15) and using equation (16) give 


m 

a 



n 


and the radiant exitance estimate at the surface is 


N 


m 
a 

n „m^ 


M(0,$) = X) E \ 

n=0 m=-n n 


( 20 ) 


Thus, the estimate of the radiant exitance depends on the coefficients of the mea- 
surement representation aJJ and the eigenvalues of the measurement operator 


IMPLEMENTATION OF TECHNIQUES 

Important to any data analysis technique is its application. Many times, ease 
of application dictates the choice of the technique. Such implementation details as 
the coordinate system, numerical techniques, conputer time and storage requirements, 
and necessary assumptions are all part of a technique and affect the final results. 


Field Representations 

The radiant exitance field at the top of the Earth's atmosphere is represented 
in this paper in two ways: in a regional grid system and in a spherical harmonic 
system. Each system has advantages. The transformations from one representation to 
the other is necessary when coitparing results. 

The first system, illustrated in figure 2, is the basic regional grid system. 

It is obtained by dividing the Earth into 10® colatitude zones and then subdividing 
each zone into an even number of regions of equal area. The equatorial zones 
80® < 0 < 90® and 90® < 0 < 100® each contain 36 regions, so that each region is 
bounded by 10® lines of colatitude and 10® lines of longitude. In all other zones, 
the number of regions is selected to make the areas of these regions as nearly as 
possible equal to those at the equator. All of these areas are referred to as either 
10® X 10® regions or simply 10® regions. 
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Three other regional grid systems are derived from this basic system. A 5® x 5° 
regional grid system is obtained dividing each 10® region in half in both colati- 
tude and longitude. A 2.5® x 2.5® regional grid system is obtained by dividing each 
5® region in half in both co latitude and longitude. A third system is obtained by 
combining adjacent 10® regions in pairs in each zone to obtain 10® x 20® regions. 

The radiant exitance field is represented in these grid systems by a piecewise con- 
stant function over the regions. 



Figure 2.- Polar view of 10® x 10® regional grid representation. 


Another way to represent the radiant exitance field is by a system of spherical 
harmonics. If X(0,<^) is a real function on the surface of a sphere, then a repre- 
sentation of X(0,^) is 


X(0,^) 



y"* (0,®) + s™ y"* (0,®) 

cn n sn 


( 21 ) 


where and are real coefficients of the spherical harmonics and 

respectively. The spherical harmonics are defined by 
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y“^(0,3') = N™ cos in® P“(cos 0) 


y“„(0,®) = n“ sin n® P”(cos 0) 




sn 


( 22 ) 


where 


n 


(2n + 1)(n - m)!(2 - 6™) 


(n + m) 1 


1/2 


Moreover/ the following properties apply: 


£ 


. dA = 

Sphere -^Sphere 




. dA = 47 c6^ 6^ 


(23) 


i 


Sphere 


y"* Y^ . dA = 0 
cn s j 


(24) 


where 6™ is the Kronecker delta function. The radiant exitance over the surface of 
the Earth can be defined hy the coefficients cJJ and sj^. 


Transformation Between Regions and Spherical Harmonics 

The regional grid system defines the radiant exitance at every point as a piece- 
wise constant function in terms of the regional values. The spherical harmonic sys- 
tem defines the radiant exitance at every point as a linear combination of spherical 
harmonics in terms of the coefficients C™ and sj^. A given set of regional values 
can be transformed to spherical harmonic coefficients hy the principle of orthogonal 
projection as follows. Multiplying equation (21) by Y^.(0,<&), integrating over the 
sphere, and taking into account the orthogonality conditions of equations (23) and 
(24) give the coefficients as 


c”" = h f X(0,®) Y‘" (0,®) dA 

Sphere 


-HI 


Since X(0,^) is known in terms of its regional values we have 
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c” = — 

n 4H 


\ I 

k=1 •'k 


y"" (0,$) dA 

cn 


Region k 


and from equation { 22 ), 


n"' K r® ^0 

C™ = ” / cos mfS / P™(cos 0) sin 0 d0 

n 471 ' k ^ n 

k=1 . 0=0, . 

k1 k1 


or 


n” k 


< ■ 3i 2 \ V\l'V 


(25) 


Similarly, 


n“ K 
k=1 


(26) 


where subscripts 1 and 2 denote the boundaries of a region. The integrals in equa- 
tions (25) and (26) are defined as 


C (^> 

m 


1'®2> = I ‘ 


cos m$ d® 


(27) 




r®2 

S (®./'^_) = / sin IDS’ d® 


(28) 


1 


1-0, .0,) ^ P”(cos 0) ,1„ 0 

n 1 2 n 


d0 


(29) 


These quantities are efficiently computed by recursive formulas (Bess et al. 1981 )• 
Equations (25) and (26) define the transformation between the regional grid system 
and the spherical harmonic system. 
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The transformation from the spherical harmonic coefficients to the regional 
values also follows from the principle of orthogonal projection. The regional values 
of a function X(0^^) are simply the average over the region, or 


f . X(0,^) dA 

_ •'Region k _ , 

X, = ^5 = A, 




Region k 


dA 


-1 j\z r®k2 


X(0,^) sin 0 d0 d^ 


(30) 


where A^ is the surface area of the kth region. Substituting for X(0,^) and 
ea^ressing the results in terms of definitions (27) through (29) yield 






This equation defines the transformation from the spherical harmonic system to the 
regional grid system. 


Relationship Between Region Size and 
Spherical Harmonic Representation 

The transformation from the spherical harmonic system to the regional system is 
accomplished by averaging over the region (eq. (30)), which smooths the representa- 
tion. Thus, the transformation from regional values back to spherical harmonic coef- 
ficients (eqs. (25) and (26)) yields a set of coefficients that represent the aver- 
aged or smoothed field. However, we wish to recapture as well as possible the origi- 
nal unsmoothed field from the regional values. Therefore, the spherical harmonic 
coefficients produced from this coordinate transformation must be increased, or 
enhanced, to represent the original field. Consider the smoothing parameter 
defined hy Pellinen (1967) and given by (appendix A) 


P 


n 


1 


J 

cos Y* 




n 


(cos Y) 


sin Y dY 


(31) 


If a field is smoothed by replacing each point value with the average value over a 
circular region with radius y*, then the coefficients and of the smoothed 

spherical harmonic representation are given by and = ^n®n* 

(1977) has shown that smoothing over square regions can be approximated with 
provided that an effective radius is used. The effective radius gives circular 
regions with the same surface area as square regions. For example, a 10° x 10® 
region corresponds to a circular region with Y* = 5.63®. The smoothing parameter 
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is presented in figure 3 for two values of Y*« Thus, to recapture the original 

n jri Jfl 

spherical harmonic coefficients and from the regional values, we first 

calculate the spherical harmonic coefficients and that represent the smooth 

field (eqs. (25) and (26)) and then divide them by that is. 


,m = e 

n '“^n ^n 


n ®n 


(32) 


This approach was followed throughout this study when spherical harmonics were com- 
puted from regional values. 



(5° X 5°) 
(10°xl0°) 


Figure 3.- Regional smoothing parauneter. 


Calculation of Important Parameters 

A radiant exitance field can be characterized in a nmber of ways. Three impor- 
tant characteristic parameters are the zonal radiant exitance, the global 
radiant exitance, and the radiant exitance gradients. the regional radiant exit- 

ance values are known, then the zonal radiant exitance Mj is 





E Ai 



AiMi 


(33) 
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where i ranges over all regions in the jth colatitude zone and where the prime 
denotes zonal values. Similarly, the global radiant exitance is 


A 

M 

g 




A 


? "A 


( 34 ) 


where A is the area of the globe and k 


ranges over all regions. It can be shown 


that the global radiant exitance is identical to Cq as given by equation (25). In 


addition, we define the pole-to-pole gradient as c!j and the equator-to-pole gra- 
dient as C 2 « The deconvolution technique produces these values directly since 
they are spherical harmonic coefficients, whereas the shape factor and parameter 
estimation techniques produce regional values from which the gradient values are 
determined hy equation (25). 


The radiation fields produced by the three analysis techniques are compared in 
two ways. The first comparison is in physical space by comparing regional values, 
and the second conparison is in spectral space by coitparing the spherical harmonic 
coefficients . 


Coitparison in physical space .- Define the difference in the ith 10® region 
between two radiation fields as The differences in the zonal radiant 

exitances follow from equation (33) with replaced by Further, we define 

the average zonal difference as 


Area weighted mean I -1 

of absolute zonal > = A ^ A*. |Am!| 
difference J ^ j=1 ^ ^ 


(35) 


We also define the average 10® regional difference as 


Area weighted mean 
of absolute 10® 
regional difference 


) 


- V’ 


k=1 


(36) 


Another descriptive parameter is 


Standard deviation of 
regional differences 


}• 


A Z A. - 


11/2 


k=1 


k k 




(37) 
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where 



E 

k=1 ^ 


We can also calculate the correlation coefficient between field 1 and field 2 on a 
10° regional scale as 


-1 




K 

E 

k=1 ^ 


1k - 


”-1 ^ ^2 

1/2 

~ _ K “ 

\^”lk “ ”l^ 


A ” E A. (M_ - M,) 
Lg k 2k 2 J 


1/2 


where 


(38) 


M 



k=1 


Comparison in spectral space ." In spectral space, we conpare the degree disper- 
sions of each field given by 


a 

n 




(n = 1, 2, N) (39) 


APPLICATION OF TECHNIQUES TO DATA SET 

The three techniques which have been presented were used to analyze wide field 
of view radiometer data to produce a radiant exitance field at the top of the atmo- 
sphere. Each technique has its own assumptions and limitations which are examined 
individually in this section. The same data set was used for all nrutierical examples, 
so that the differences in the derived radiation fields are entirely due to the tech- 
nique and the assumptions employed. 


Earth Radiation Data 

The measurement data used in this study to define the Earth's longwave radiation 
field were obtained from the Earth Radiation Budget (ERB) instrument aboard the 
Nimbus 6 satellite. A description of the ERB instrument and calibration is given by 
Smith et al. (1977). The data tapes were supplied by the National Oceanic and Atmos- 
pheric Administration. 
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The ERB instrument obtained both fixed wide-angle and scanning narrow-angle 
measurements. The data considered here are tbe fixed wide-angle measurements at 
satellite altitude. One data channel recorded the total Earth measurement (0.2 to 
50 ^lm) and another channel recorded the shortwave measurement (0.2 to 3.8 Mm). The 
longwave contribution is the difference between these two measurements. The ERB 
instrument operated with a duty cycle of 2 days on and 2 days off. Measurements were 
taken every 16 seconds. The data set for this study consists of approximately 46 000 
wide-angle measurements taken during August 1975. The emphasis in this stu<^ is not 
on the measurement data or the derived fields but on the differences in the fields 
derived with different assumptions and techniques. 

The model of the radiant exitance field is independent of time; thus all mea- 
surement data were collected into 5® x 5® regions at satellite altitude (h = 1070 km) 
and averaged. We processed only 1664 averages instead of the 46 000 individual 
measurements . 


Shape Factor Technic[ue 

With the shape factor technique, the radiant exitance at the top of the atmo- 
sphere is derived from each measurement by dividing by a shape factor (eq. (6)). The 
two basic assumptions are that ( 1 ) the radiant exitance at the top of the atmosphere 
is constant over the field of view and (2 ) the estimate of radiant exitance is asso- 
ciated with the region which contains the nadir point. As a result of the first 
assumption, which is a great simplification, the technique yields point estimates of 
the radiant exitance at nadir that are about 5 percent in error (Weaver and Green 
1980, Green 1981). However, these estimates are seldom interpreted as point esti- 
mates but are used to represent the average radiant exitance over an area. In fact, 
these estimates best represent the average radiant exitance over a circular surface 
area centered at nadir with a 20® diameter. On the average with the shape factor 
technique, the average radiant exitance over this area can be estimated to within 
about 1 percent (Weaver and Green 1980, Green 1981). In this study we do not con- 
sider point estimates or estimates over circular regions, but estimates over square 
regions, over zones, and over the globe. All estimates of radiant exitance whose 
nadir points are in a particular area are averaged. Thus, the assumption of constant 
radiant exitance over the field of view is not as restrictive as one might imagine. 
The larger the area, the more appropriate the assumption. The second assiomption of 
associating the estimate of radiant exitance with the region that contains the nadir 
point follows from the geometry of the measurement. For a uniform field, an incre- 
mental area centered at nadir has more influence on the measurement than any other 
incremental area in the field of view. 


In addition to these two assumptions, two simplifying assximptions have been made 
about the directional function. First, we assume that R(9) is the same for all 
points on the globe. This assumption was shown to yield a shape factor which is 
independent of the directional function, namely, F = [R^/(Rg + h)]^. Second, we 
model h as a constant, 1070 km for the Nimbus 6 data set. Thus, for an Earth 
radius of 6378 km and an atmospheric height of 30 km, we have a shape factor of 


\6378 + 30 + 1070 


f- 


0.7343 


(40) 
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With these assumptions the only parameter that affects the derived radiation 
field is the size of the regional grid system that is employed. The 46 000 measure- 
ments have been averaged over a 5® x 50 gj^xd system at satellite altitude. Each of 
these averages is divided by the shape factor to produce an estimate of the average 
radiant exitance of a 5® 5® region at the top of the atmosphere. These estimates 

can easily be combined to yield 10 ® ^ 10 ® regional radiant exitance estimates. 

Another approach would be to average the 46 000 measurements over a 10® x 10® grid 
system at satellite altitude and then to reduce each average to the top of the atmos- 
phere. This approach allows for less spatial variation and can be considered infe- 
rior to the first approach. The difference between these two solutions illustrates 
the effect of the grid system size. 


Table I presents these two solutions and the differences between them. The 
global radiant exitances as given by eqpiation (34) differ by 0.30 W-m"^ because of 


TABLE I.- SOLUTIONS FROM SHAPE FACTOR TECHNIQUE ILLUSTRATING 
EFFECT OF GRID SYSTEM 



Case 1, 

Case 2, 

Case 1 

Parameter 

data averaged 

data averaged 

minus 


over 5® regions 

over 10® regions 

case 2 

-.9 

Global radiant exitance, W-m 

235.39 

235.09 

0.30 

Pole-to-pole gradient, W-m 

11.40 

11.61 

-.21 

—2 

Equator-to-pole gradient, W-m 

-22.03 

-22.27 

.24 

“2 

10® regional radiant exitance, W-m : 




Pacific, tropical 

266.03 

265.58 

.45 

Pacific, high latitude 

191.33 

191.23 

.10 

Atlantic, subtropical 

268.29 

268.02 

.27 

Sahara 

279.59 

281.60 

-2.01 

Greenland 

221.09 

220.71 

.38 

South Pole 

124.50 

118.57 

5.93 

_2 

Area weighted mean of absolute zonal differences, W-m 


0.34 

—2 

Area weighted mean of absolute 10® regional differences, W-m 

.83 

_2 

Standard deviation of 10® regional differences, W-m ... 


1.44 


the grid system. The pole-to-pole gradient and the equator-to-pole gradient are 
given by the spherical harmonic coefficients and respectively, at the top 

of the atmosphere, as determined from equations (25) and (32). The global value Cq 
can also be calculated from these equations which give the same value as equa- 
tion (34), since the smoothing is given by Pg ” ^ • Also included in table I are the 
10® regional radiant exitances for six specific regions. The average zonal differ- 
ence as computed from equation (35) is 0.34 W-m""^, nearly the same as the global 
difference of 0.30 w-m" . In general, the smaller the area in question, the greater 
the difference. This is illustrated by the average regional difference of 0.83 W-m"^ 
as confuted from equation (36) and their standard deviation of 1.44 W-m"^ from equa- 
tion (37). If we consider the range of regional differences to be ±2Cf about their 
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—2 —2 

mean/ then they range over 0,30 ± 2 (1«44), or from -2.58 W-m to 3.18 W-m • (The 
mean regional difference is given by the global difference (0.30 w-m“^).) The six 
specific regional differences all fall within this range of differences except for 
the South Pole. 

The zonal differences for the two solutions are plotted in figure 4. Notice the 
unusually large difference for 160® < 0 < 170® compared with the small differences 
for all other zones. The problem is due to nonuniform satellite sampling combined 
with a large latitudinal gradient. Using the 5® grid system, we can determine the 
5® zonal radiant exitance from 160® < 0 < 165® to be 145.38 W-m^ and from 
165® < 0 < 170® to be 123.35 W-m”^. By averaging the two 5® zonal radiant exitances 
with the appropriate area weighting, we obtain 136.16 W-m which corresponds to the 
5® grid system solution. In the 10® grid system, this 5® zonal distinction is not 
made; averaging all data between 160® < 0 < 170® results in 129.54 W-m“^. If the 



Figure 4.- Zonal differences due to grid system (5® or 10®) for 

shape factor technique . 


satellite sample were uniform with area, then the 10® zonal results would be the 
same. Since the number of measurements for 160® < 0 < 165® was 841 and the number of 
measurements for 165® < 0 < 170® was 2154, the southern zone was weighted more. This 
nonuniform sampling is a result of orbital geometry and data elimination due to Sun 
contamination (Bess et al. 1981). The difference in the 160® < 0 < 170® zonal rad- 
iant exitance is also a function of the latitudinal gradient. If the 5® zonal values 
were equal, then the number of measurements would not matter. Notice that the dif- 
ference in the 10® < 0 < 20® zonal radiant exitance is not nearly as large. This 
results from a smaller latitudinal gradient for the northern zone. Nonuniform sampl- 
ing combined with a large latitudinal gradient is also the cause of the large differ- 
ence in radiant exitance for the South Pole 10® region. 

The question arises as to why large differences in radiant exitance do not occur 
for the 170® < 0 < 180® zone. The Nimbus 6 orbit has an inclination slightly less 
than 100®, which means that the nadir point slightly exceeds a latitude of ±80®. 

Thus, we have sampling for the 170® < 0 < 180® zone, although all of it is near 

0 = 170®. A similar situation exists at the North Pole. For the 5® grid system, the 

satellite samples 170® < 0 < 175® but not 175® < 0 < 180®. Thus, we must extrapolate 

into the 175® < 0 < 180® zone. The rule has been adopted to set this zonal radiant 
exitance equal to the 170® < 0 < 175® estimate. Since both 5® zones have the same 
radiant exitance, we understandably obtain an estimate for the 5® grid system which 
is similar to the 10° grid system estimate. 
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If reducing the size of the grid system from a 10® grid to a 5® grid changes the 
160® < 0 < 170® zonal estimate, then what effect does further reducing the grid sys- 
tem to 2,5® have? The 5® grid, like the 10® grid, overweights the estimates toward 
the polar values because of the satellite sampling. The 2.5® grid reduces this bias 
and increases the zonal radiant exitance to 137.92 W-m“^ for 160® < 0 < 170®. How- 
ever, this small change may not be significant conpared with the accuracy of the 
shape factor technique. 


Parameter Estimation Technique 

With the parameter estimation technique, all measurements are processed together 
as a batch and the radiant exitance estimate at the top of the atmosphere is a least 
squares fit to the data using equation (14). The two basic assumptions are that 

( 1 ) the radiant exitance at the top of the atmosphere is constant over a region and 

(2) the directional function is known for each region. The first assumption is rea- 
sonable since we desire to estimate the mean radiant exitance over a region, and the 
estimates we obtain by assuming the radiant exitance to be constant closely resemble 
the mean regional radiant exitances. The second assumption of known directional 
functions introduces errors into the estimates. Estimates of the radiant exitance 
for large areas like the globe or a zone are not as sensitive to the directional 
functions as regional estimates. These errors are examined in this section 1:^ a 
numerical example. 

One problem with implementing the parameter estimation technique (eq. (14)) is 
the size of the matrix. There are 416, 10® 10® regions over the globe. 

Thus, the 416 x 416 B*^B matrix contains 173 056 elements or 83 736 if symmetry is 
taken into account. The required conputer storage was reduced hy one-fourth by esti- 
mating the radiant exitance in 10® ^ 20® regions. These 10® x 20® estimates are 

divided into 10® ^ io® estimates by a curve fitting scheme presented in appendix B. 

A source of error is the choice of the directional function for the 208 rectan- 
gular 10® X 20® regions. This function has been denoted is included in 

the calculation of B^^^ (eq. (12)), The effect of different functions was excimined 
by two numerical examples in which all 10® ^ 20® regions were assumed to have the 
same directional function. In the first case the radiation was assumed to be 
Lambertian, or Rj^(Q) = 1. In the second case a directional function with a con- 
siderable amount of limb darkening was assumed and is given by 

R(0) = 1.106 035)^0.0547(1 - sec 9) + 0.09375 {exp [2 .24 (1 - sec 9)] - 1 (41) 

where 0® < 0 < 90®. This function represents a lower bound to a family of limb dark- 

ening functions based on Nimbus 2 data (Raschke et al. 1973). If these two direc- 
tional functions are considered upper and lower bounds, the two solutions establish 
bounds on the radiant exitance due to uncertainty in the regional directional func- 
tions. 


Part of the implementation of the parameter estimation technique is the selec- 
tion of a quadrature formula of the influence coefficients given by equation (12), 
that is. 
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where is the influence of the kth region (10® x 20® region) on the ith measure- 

ment/ and the integration is over the ith field of view that is also in the kth 
region. A very practical numerical integration scheme is to subdivide the regional 
grid system into smaller siabregions and approximate the integral of equation (12) by 
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where j ranges over all the subregions that are in both the ith field of view and 
the kth region. The subscript j denotes that a parameter is evaluated at the 
center point of the jth subregion. Also, a siibregion is included in the summation if 
its center point qualifies and is excluded otherwise. This means that the boundaries 
of the integration are approximated 1:^ meridional and latitudinal lines. This bound- 
ary error and the quadrature error are reduced by making the subgrid system smaller 
at the expense of computational effort. Numerical examples are given for two differ- 
ent subgrid systems. 

Three numerical examples are presented in table II. The measurement data for 
each example are the same and were averaged into 5® x 5® regions at satellite alti- 
tude. This reduced the number of measurements to 1664 values, which were processed 


TABLE II.- SOLUTIONS FROM PARAMETER ESTIMATION TECHNIQUE ILLUSTRATING 
EFFECTS OF QUADRATURE SUBGRID SYSTEM AND DIRECTIONAL FUNCTION 


Parameter 

Case 1, 

Lambertian R, 
5® quadrature 

Case 2, 

Lambeirtian R, 
2.5® quadrature 

Case 3, 

limb darkening R, 
2.5® quadrature 

Case 1 
minus 
case 2 

Case 2 
minus 
case 3 

Global radiant exltance, W-m 

235.37 

235.42 

235.42 

-0.05 

0 

Pole-to-pole gradient, W-m 

11.59 

11.61 

11.58 

-.02 

.03 

_2 

Equator-to-pole gradient, W-m 

-23.44 

-23.35 

-23.21 

-.09 

-.14 

«2 

10* regional radiant exitance, W-m : 






Pacific, tropical 

270.72 

270.64 

270.73 

.08 

-.09 

Pacific, hi^ latitude 

187.63 

186.87 

187.72 

.76 

-.85 

Atlantic, subtropical 

279.65 

280.04 

278.91 

-.39 

1.13 

Sahara 

307.34 

308.38 

305.45 

-1.04 

2.93 

Greenland 

225.30 

226.29 

225.59 

-.99 

.70 

South Pole 

97.98 

96.55 

99.71 

1.43 

-3.16 

Con^uter time , sec 

120 

335 

456 



Area vreighted mean of absolute zonal differences, W-m”^ • 



0.56 

0.85 

Area weighted meem of absolute 10» regional differences. 

W-m“^ 


.76 

1.20 

Standard deviation of 10® regional differences, W-m~^ •• 



1.09 

1.51 
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with the parameter estimation technique (eq. (14)). In the first two cases / we 
assumed the radiation to be Lambertian and used a 5® and a 2.5® subgrid quadrature 
scheme. The average zonal difference is 0.56 W-m and the average regional differ- 
ence is 0.76 w-m"^. Since the global difference is approximately zero (-0.05 W-m"^)^ 
the average zonal and regional differences reflect the true differences in the shape 
of the radiant exitance field. It is reasonable to ejq)ect the 2.5® subgrid q[uadra- 
ture scheme to give the better results since the quadrature errors are smaller. The 
computer time^ also listed in the table, is about three times greater for the smaller 
subgrid quadrature. However, 335 seconds of coitputer time on a Control Data CYBER 
173 computer is not considered excessive, so that a 2.5® quadrature scheme was 
adopted for this study. The zonal differences are plotted in figure 5 and are 
relatively small except for the northern and southern zones. 



Figure 



Cases 2 and 3 in table II represent the differences due to different directional 
functions. The differences are very small on a global scale, larger on a zonal 
scale, and still larger on a regional scale. The global radiant exitance was the 
same for both cases, and the average zonal difference was 0.85 W-m"*^. The average 
regional difference was 1.20 W-m and we would expect the individual regional dif- 
ferences to range between ±3.02 W-m”^. Also, assuming a limb darkening function 
increases the computer time from 335 to 456 seconds. It is not clear which of these 
two solutions best represents the true radiant exitance field. The first assumes no 
limb darkening, while the second probably overestimates the limb darkening. Assximing 
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a limb darkening function between these two would probably give better results. 
Moreover, if the directional character of each region were known then the solution 
would be improved even further. The purpose here, however, is to define on a global, 
zonal, and regional scale the differences due to the extreme directional functions. 


Deconvolution Technique 

With the deconvolution technique, we take advantage of the fact that spherical 
harmonics are the eigenfunctions of the measurement operator. If the measurements 
are used to define the radiant exitance field at satellite altitude in terms of a 
spherical harmonic expansion, then the field is easily reduced to the top of the 
atmosphere by the eigenvalues (eq. (20)). This solution is based on the following 
three simplifying assumptions: (1) the sensor integrates incoming radiation over its 

field of view with its directional response S(a) being a function only of the nadir 
angle, (2 ) the directional dependence of exiting radiation at each point of the top 
of the atmosphere is a function only of the zenith angle of the exiting ray, and 
(3) continuous measurements are available over a sphere concentric with the Earth. 

The first assumption is applicable to most wide field of view sensors because of 
their design. Throughout this study we have used S(a) = cos oc, or a flat-plate sen- 
sor model. The second assumption, which states that R(©,^, = R(0), is a reason- 
able approximation for Earth emitted radiation. The third assumption contains two 
requirements. The first is that the measurements be continuous over a sphere. Since 
physical systems provide discrete measurements, we calculate the spherical harmonic 
coefficients with equations (25) and (26) and then approximate the continuous case 
with the smoothing parameter The second requirement is that the orbit be 

circular. 

The implementation of the deconvolution technique requires that we define the 
directional function R(^) and that we truncate the spherical harmonic expansion at 
some finite degree. To investigate the dependence of the solution on the directional 
function, the limiting cases of a Lambertian function (R(Q) = 1) and a limb darkening 
function (eq. (41)) have been examined. These directional functions affect the solu- 
tion through the eigenvalues (eq. (17)), which are given in table III for degree 0 
through degree 12, Note that = 0.7343 for both directional functions; thus, the 


TABLE III.- EIGENVALUES OF MEASUREMENT OPERATOR 
[R^ * 6378 + 30 = 6408 km; h = 1070 km] 


Degree, 

n 


for - 

Lambertian R 

Limb darkening R 

0 

0.7343 

0.7343 

1 

.7217 

.7233 

2 

.6975 

.7019 

3 

.6632 

.6713 

4 

.6208 

• 6331 

5 

.5726 

.5893 

6 

.5214 

.5419 

7 

• 4693 

• 4929 

8 

• 4185 

.4443 

9 

.3707 

.3974 

10 

• 3267 

• 3535 

11 

.2874 

• 3132 

12 

.2526 

.2770 


23 




global radiant exitance is unaffected by directional function. It is also interest- 
ing that is the shape factor (eg. (40)) used in the shape factor technique. To 

examine the effect of the degree of truncation on the solution, we have expanded the 
measurements at satellite altitude in a series of spherical harmonics to degree 24 by 
ecjuations (25) and (26) and approximated the continuous solution by using based 

on 5® ^ 5® regions. The degree dispersion (eq. (39)) for this solution is given in 
figure 6. The radiant exitance estimate at the top of the atmosphere was confuted 



Figure 6.- Degree dispersion of deconvolution solution based on 
Lambertian directional function. 


by dividing the spherical harmonic coefficients by the eigenvalues based on 

Lambertian radiation. The divergence of these two solutions in figure 6 is due to 
the smoothing effect of the measurement operator. It is seen that the radiant 
exitance spectrum measured at satellite altitude can be divided into two regions, one 
in which log decreases approximately linearly with n for n < 15, and one in 

which log is rather constant for n > 15. The spectrum at the top of the atmo- 
sphere likewise consists of two regions, one in which log decreases somewhat 
linearly for n < 15, and one in which log increases somewhat linearly for 
n > 15. The nearly constant spectrum at satellite altitude is interpreted as mea- 
surement noise and produces the divergence of the solution at the top of the atmo- 
sphere since the eigenvalues are a decreasing sequence of values. Thus, it is 
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concluded from figure 6 that the coefficients for n > 15 are due primarily to noise 
and contain little information concerning the radiant exitance field; the expansion 
is therefore truncated at N = 15, 

The interpretation of a nearly constant spectrum as random noise is illustrated 
by examining the spectrum of a random field. This field is defined on a 5® regional 
grid system where each regional value is a normal random variable with mean of zero 
and standard deviation of 10 W-m""^. The expected value of the spectrum for this 
field is derived in appendix C and presented in figure 6. Although this spectrum for 
a random field gradually increases for n > 15, it gives support to the interpreta- 
tion of a nearly constant spectrxim as noise. 

To investigate the sensitivity of the solution to the directional functions and 
the degree of truncation, three cases were examined, as given in table IV. Two 
limiting directional functions were examined and two different truncation degrees. 


TABLE XV.- SOLUTIONS FROM DECONVOLUTION TECHNIQUE ILLUSTRATING DEGREE OF 
TRUNCATION AND OF DIRECTIONAL FUNCTION 



Case 1 , 

Case 2, 

Case 3, 

Case 1 

Case 2 

Parameter 

Lambertian R, 

Lambertian R, 

limb darkening R, 

minus 

minus 


degree 15 

degree 12 

degree 12 

case 2 

case 3 

Global radiant exitance, W-m 

235.39 

235,39 

235.39 

0 

0 

Pole-to-pole gradient, W-m 

11.60 

11.60 

11.58 

0 

.02 

Equator-to-pole gradient, W-m“^ 

-23.19 

-23.19 

-23.05 

0 

-.14 

10° regional radiant exitance, W-m”^: 






Pacific, tropical 

259.70 

263.94 

264.26 

-4.24 

-.32 

Pacific, high latitude 

190.99 

191.40 

191.65 

-.41 

-.25 

Atlantic, subtropical 

281.28 

274.47 

273.52 

6.81 

.95 

Sahara 

283.54 

289.95 

288.46 

-6.41 

1.49 

Greenland 

215.67 

213.77 

214.32 

1.90 

-.55 

South Pole 

1 109.38 

111.05 

113.15 

-1.67 

-2.10 

_2 

Area weighted mean of absolute zonal differences, W-m 



0.85 

0.55 

Area weighted meein of absolute 10® regional differences 

, W-m“2 


3.94 

.75 


Standard deviation of 10® regional differences, W-m” 


4.89 


.94 


As discussed previously, a reasonable choice of the degree of truncation is N = 15 
based on the present 1-month data set. However, Bess et al. (1981) have investigated 
12 months of Nimbus 6 data and concluded that N = 12 is the appropriate e3q>ansion 
limit. Therefore, the two solutions based on N = 12 and 15 are compared in 
table IV. This comparison reveals that the average regional difference due to trun- 
cation is 3.94 W-m" • Also, the average regional difference due to different direc- 
tional functions is 0.75 W-m" . Thus, we conclude that on a regional scale, the 
degree of truncation has more influence on the solution than the directional func- 
tion. The effect of truncation is reduced, however, on a zonal scale, since the 
addition of the high-frequency coefficients (n = 13, 14, 15) does not affect the 
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zonal values as much as the regional values* The average zonal difference due to 
truncation is 0.85 and due to different directional functions is 0.55 W-m" • 

The individual zonal differences are plotted in figure 7. The degree at which the 
spherical harmonic expansion should be truncated is very important to the analysis. 
However, there is not a clear criterion for choosing this degree. For the remainder 
of this study we use a spherical harmonic expansion to degree 12 for the deconvolu- 
tion technique. 
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(a) Difference due to truncation 
(case 1 minus case 2). 

\ I L. 



Figure 7.- Zonal differences for deconvolution technique. 


COMPARISON OF TECHNIQUES 

The three data analysis techniques have been derived and examined individually 
to establish their necessary assumptions and characteristics. In this section, these 
techniques are compared numerically. Since each technique processes the same mea- 
surement data, any differences in the derived radiation fields are a result of the 
techniques themselves • 

A summary of the three radiation fields that result from the three techniques is 
presented in table V and the zonal differences are plotted in figure 8. The shape 
factor solution corresponds to the 5® regional grid system and was presented in 
table I. This case was chosen to represent the shape factor technique so that all 
techniques would have not only a common data set but also a common grid system. The 
parameter estimation solution corresponds to the Lambertian directional function and 
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TABLE V.- SOLUTIONS FROM THREE TECHNIQUES FOR COMPARISON 


Parameter 

Shape 

factor 

Parameter 

estimation 

Deconvolution 

Shape factor 
minus 

parameter estimation 

Shape factor 
minus 

deconvolution 

Parameter estimation 
minus 

deconvolution 

Global radiant exitance r W-m"^ 

235.39 

235.42 

235.39 

-0.03 

0 

0.03 

Pole-to-pole gradient, 

11.40 

11.61 

11.60 

-.21 

-.20 

.01 

Equator-to-pole gradient, W-m“^ 

“22.03 

“23.35 

“23.19 

1.32 

1.16 

-.16 

10* regional radiant exitance# W-m : 







Pacific, tropical 

266.03 

270.64 

263.94 

-4.61 

2.09 

6.70 

Pacific, high latitude 

191.33 

186.87 

191.40 

4.46 

-.07 

“4.53 

Atlantic, subtropical 

268.29 

280.04 

274.47 

-11.75 

-6.18 

5.57 

Sahara 

279.59 

308.38 

289.95 

-28.79 

-10.36 

18.43 

Greenland 

221.09 

226.29 

213.77 

-5.20 

7.32 

12.52 

South Pole 

124.50 

96.55 

111.05 

27.95 

13.45 

“14.50 

—2 

Area weighted mean of absolute zonal differences, W-m ... 


7.42 

4.04 

4.84 

Area weighted mean of absolute 10* regional differences, W- 



10.65 

5.95 

8.60 

Standard deviation of 10® regional differences. 

W-m"2 


13.51 

7.44 

10.78 

Correlation (regional scale) 




.95 

.98 

.97 









a 2.5® quadrature scheme and was first presented in table II. The Lambertian func- 
tion was chosen for its simplicity, and the 2.5® quadrature scheme for its accuracy. 
The deconvolution solution taken from table IV corresponds to a spherical harmonic 
expansion to degree 12 and a Lambertian directional function. The 12th degree expan- 
sion was chosen since it probably represents the appropriate e3q>ansion limit for 
monthly information over a yearly cycle, as discussed by Bess et al. (1981). The 
Lambertian function was chosen so that the parameter estimation and the deconvolution 
would have no differences due to the directional function. All three techniques 
produced essentially the same estimate of the global radiant exitance. This is char- 
acteristic of a flat -plate sensor and the wide field of view measurement. The other 
pertinent paraimeters, however, depend on the technique employed. 

The differences in these three solutions also are presented in table V. The two 
solutions that correlate best on a 10® regional scale (eq. (38)) are the shape factor 
and the deconvolution solutions, their correlation being 0.98. Their average zonal 
difference is 4.04 W-m and their average regional difference is 5.95 w-m • It is 
surprising that the two techniques with the highest correlation should vary by this 
large amount. The regional standard deviation of 7.44 W-m*"^ implies that about 
95 percent of the regional differences are between ±15 W-m . The two solutions with 
the least correlation are the shape factor and the parameter estimation solutions. 
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SHAPE FACTOR MINUS PARAMETER ESTIMATION 


PARAMETER ESTIMATION MINUS DECONVOLUTION 



SHAPE FACTOR MINUS DECONVOLUTION 



COIATITUDE, deg 


Figure 8,- Zonal differences for three techniques. 


_2 

The regional standard deviation for these two solutions was 13.51 w-m . The distri- 
butions of these differences are presented in histograms in figure 9. These distri- 
butions are approximately normal and centered about zero. The common global value 


NUMBER 

OF 

REGIONS 



REGIONAL DIFFERENCE. W-m“^ 


Figure 9.- Histograms of regional differences for the three techniques* 


for the three solutions results in the zero means. Again, we see that the shape 
factor and the deconvolution solutions are most similar on a regional scale. 

The three solutions can also be compared hy examining their spatial spectra. 

The basic representation of the deconvolution solution is in the spherical harmonic 
coordinate system. The shape factor and the parameter estimation solutions have been 
transformed to the spherical harmonic system, and the spatial spectrums of all three 
solutions are conpared in figure 10. The spectrum of the parameter estimation solu- 
tion decreases up to about degree 15 and then exhibits a different type of behavior. 
The deconvolution spectrum follows the parameter estimation spectrum very closely to 
degree 12 and is zero from there on because of the 12th degree truncation. The 
shape factor spectrum gradually departs from the other two spectrums, the departure 
becoming significant after degree 3. 
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Figure 10.- Degree dispersion for the three techniques. 


By comparing the spectral content of the three solutions, we conclude that the 
parameter estimation and deconvolution solutions are most similar on this basis to 
degree 12 and that the shape factor solution does not compare favorably with the 
other two solutions. Recall that the deconvolution technique enhances the high- 
frequency spatial components of the solution by dividing by a decreasing sequence of 
eigenvalues (eq. (20)). The shape factor technique, on the other hand, divides all 
measurements by a single shape factor (eq. (6)) which is equivalent to dividing all 
frequency components by the first eigenvalue \q . Thus, the difference between the 
two spectra increases with degree. Since the wide field of view sensor smooths the 
measurements, it is reasoned that the deconvolution technic[ue yields the more accu- 
rate solution. Although not e^licit in the formulation, it can be concluded that 
the parameter estimation solution also enhances the high-frequency components, since 
it is so similar to the deconvolution solution. This enhancement was also demon- 
strated by Smith et al. (1975). Moreover, the 10® x 20® grid system does not seem to 
degrade the parameter estimation solution. It is remarkable that the spectra of the 
parameter estimation and deconvolution solutions are so similar considering that the 
spherical harmonic coefficients of the deconvolution solution were initially computed 
at satellite altitude with 5® x 5® regional data and the coefficients of the param- 
eter estimation solution were coir 5 >uted at the top of the atmosphere with 10° x 20® 
regional values. 
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However, the parameter estimation and the deconvolution solutions do not coit 5 )are 
favorably^on a regional scale (from table V, the average regional difference is 
8.60 w-m ^ ) • The problem results from the coordinate system transformations . The 
deconvolution solution is defined in terms of spherical harmonics. The parameter 
estimation solution is defined on a 10® ^ 20® grid system, converted to a 10® x 10® 
grid system, and then converted to a spherical harmonic representation.. Tradition- 
ally, the transformation from regions to spherical harmonics has been considered 
valid to a degree number N equal to 180 divided by the region size, or N = 18 for 
10® regions. Rapp (1977), however, has shown that this transformation can be 
expanded even beyond this degree. We have extended the transformation to N = 24 in 
figure 10 but have only compared to degree 12, since that is the limit of the decon- 
volution spectrum. Certainly the transformation from regions to spherical harmonics 
yields a spectrum of the parameter estimation solution reliable to degree 12. This 
is further substantiated by the favorable comparison of the two spectra. The 
regional comparison, on the other hand, necessitated transforming the deconvolution 
solution from spherical harmonics to regional values. Traditionally, this transfor- 
mation is valid for grid sizes equal to 180 divided by N, or 15® regions for the 
12th degree spherical harmonic e:q>ansion. Thus, the 12th degree spherical harmonic 
solution is not adequate to establish 10® regional values. Traditionally, a spher- 
ical harmonic e3q>ansion to degree 18 would be required to yield 10® regional values. 
Rapp finds that even an 18th degree expansion is inadequate and suggests that a 27th 
degree spherical harmonic expansion would be necessary to recover 10® regional 
values. Thus, the present deconvolution solution of degree 12 is quite insufficient 
to yield 10® regional values. For this reason the parameter estimation and decon- 
volution solutions do not compare favorably on a 10® x io® grid system. Since only a 
12th degree solution is available to define the regional values, they are smoothed 
because of the lack of the high-frequency terms. As previously discussed, the shape 
factor solution is also a smoothed representation, since it does not enhance the 
high-frequency components of the solution. Thus, these two smooth solutions, from 
the shape factor and deconvolution techniques, compare favorably on a regional basis. 
Nevertheless, the parameter estimation technique may give more reliable 10® x iq® 
regional values than either of these two smoothed solutions. However, because the 
true. field is unknown, one cannot demonstrate the accuracy of any analysis technique 
on the basis of these results. 

As previously discussed, use of 10® x 20® regions in the parameter estimation 
technique did not seem to degrade the spectrum of the solution, at least not to 
degree 12. From the deconvolution technique we can obtain the grid size resolution 
by the traditional method of dividing 180® by the degree number 12. This would imply 
that 15® X 15® regions, of which there are 185 to cover the globe, can be resolved. 
However, there are 208 regions of size 10® x 20® over the globe. Thus, we conclude 
that a 10® x 20® regional system is adequate. 

The three numerical excimples considered in this section are summarized in 
table VI. The sensitivity of the three techniques to the assumed directional func- 
tion represents the difference in the solution using the Lambertian and limb darken- 
ing functions. The parameter estimation technique is more sensitive than the decon- 
volution technique. The shape factor technique is independent of directional 
function. The computer storage requirements and running times are also given for 
conparison and should be interpreted relatively, since the skill of the programmer 
and the type of computer affect these numbers. We can conclude, however, that the 
coirputational burden for the shape factor and deconvolution techniques is small 
compared with that for the parameter estimation technique. Since the quadrature 
scheme and the directional function greatly affect the coitputational burden, more 
detailed entries are given in table VI for the parameter estimation technique. 
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TABLE VI.- SUMMARY OF NUMERICAL EXAMPLE SOLUTIONS 


Technique 

Shape factor 

Parameter estimation 

Deconvolution 

Algorithm 


M = 

Hte.9) = Z h 




n *=0 m=-n 

Data processing 

Individual 

Batch 

Batch 

Data set 

Averaged over 5® x 5® 

Averaged over 5® x 50 

Averaged over 5® x 50 regions 


regions 

regions 


Grid system 

5® X 50 regions 

10 ® X 20 ® regions 

Spherical harmonics 

Directional function, R 

Homogeneous 

Lambertian or limb darkening 

Lambertian or limb darkening 

AssuJJ^Jtions 

1- Radiant exitance constant 

1 • Radiant exitance constant 

1. Sensor response function 


over FOV 

over 5® x 5 ® regions 

of a only 


2. Estimate associated with 

2. R known for each grid 

2. R function of 9 only and 


nadir point 

region 

homogeneous over globe 



3. Extrapolate 10® x 20 ® 

3. Data defined on sphere 



region to 10 ® x 10 ® 
region 

4. Solution truncated at N = 12 

Sensitivity to 




directional function,^ W-m 




Global difference 

None 

None 

None 

Zonal difference 

None 

0.85 

0.55 

10 ® regional difference 

None 

1.20 

0.75 

Computer storage 

17 000 

^92 000 ^72 000 

21 000 

Computer running 

1 

L D; ^456 *^151 

5 

time, sec 


Lamb; ^335 ^^120 



^Difference between solutions with Lambertian and with limb darkening directional functions. 
“For quadrature scheme with 2.5® subgrid. 

^For quadrature scheme with 5® subgrid. 


Finally, the advantages and disadvantages of the three techniques are summarized 
in table VII. The shape factor technique requires a minimim of computational time, 
is independent of the directional function, and does not require knowledge of the 
complete data set. The deconvolution technique is fast computationally, gives accu- 
rate spherical harmonic coefficients, and produces no quadrature errors. They both, 
however, have the disadvantage of inaccurate regional values. Other than the conpu- 
tational burden, the parameter estimation technique is impressive, yielding accurate 
spherical harmonic coefficients from 10® x 10® regional values. This implies that 
the regional values are reliable. Because the directional function can be defined at 
each quadrature subregion, the technique has considerable flexibility. 


table VII.- ADVANTAGES AND DISADVANTAGES OF TECHNIQUES 


Technique 

Advantages 

Disadvantages 

Shape factor 

Simple, fast 

Requires little storage 

Popular 

Does not require complete data set 

Accurate global mean 

Independent of R 

Responds to local discontinuity 

Does not enhance high frequency 
Inaccurate regional values 
Inaccurate high-frequency 

spherical harmonic components 

Parameter estimation 

Enhances high frequency 
Heterogeneous R over globe 
Accurate spherical harmonic 
coefficients 

Requires large storage 
Requires conplete data set 
Extrapolation from 10® x 20® 
regions to 10 ® x 10 ® regions 
Quadrature errors 

Deconvolution 

Fast 

Accurate spherical harmpnic 
coefficients to degree 12 
Requires little storage 
Enhances high frequency 
No quadrature errors 

Inaccurate regional values 
Requires conplete data set 
defined on sphere 
R function of 0 only 
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CONCLUSIONS 


Three data analysis techniques have been examined for estimating radiant exit- 
ance at the top of the atmosphere from Earth radiation budget observations. The 
shape factor r parameter estimation, and deconvolution techniques have been compared 
by examining their assiimptions and characteristics, applying them to the scime data 
set, and transforming their results to the same coordinate system. 

The characteristics of the shape factor technique were examined by estimating 
the radiation field over a 5® and a 10® grid system. The difference between these 
two approaches was 0.30 W-m~^ on a global scale and 0.83 W-m"^ on a 10® regional 
scale. These results are independent of the directional function since we assumed 
the same directional function for all points on the globe. This assxamption was shown 
to yield a shape factor which is independent of the directional function. 


Two factors affect the radiation field derived by the parameter estimation tech- 
nique: the directional function and the quadrature formula. The difference between 

using a Lambertian and a limb darkening function was negligible on a global scale and 


—2 

1.20 W-m on a 10® regional scale, 
ture scheme was 0.05 w-m ^ 


The difference between a 2.5® and a 5® quadra- 
on a global scale and 0.76 W-m ^ on a 10® regional scale. 


Implementation of the deconvolution technique requires that we define the direc- 
tional function and that we truncate the spherical harmonic e>q>ansion at some finite 
degree. Neither one of these two factors affected the global average. On a 10® 
regional scale, the difference between using a Lambertian and a limb darkening func- 
tion was 0.75 w-m” and the difference between truncating the spherical harmonic 
expansion at degree 12 and at degree 15 was 3.94 W-m Thus, on a regional scale, 

the degree of truncation has more influence on the derived radiation field than the 
directional function. 


The three data analysis techniques have been inter compared. All three tech- 
niques produced essentially the same estimate of the global radiant exitance. The 
two solutions that have the highest correlation on a 10® regional scale are the^shape 
factor and deconvolution solutions. Their 10® regional difference is 5.95 w-m 
The favorable comparison of these two solutions is due to their smoothness. The 
shape factor smooths the solution by not enhancing the high spatial frequencies. The 
deconvolution technique smooths the solution by truncating the spherical harmonic 
representation at a finite degree. 

Another comparison of the three solutions is made by examining their spatial 
spectra. The parameter estimation and the deconvolution spectra are very close out 
to degree 12 where the deconvolution spectrum is truncated. The shape factor spec- 
trum gradually departs from these two spectra, the departure becoming large after 
degree 3. The parameter estimation and the deconvolution solutions do not compare 
favorably on a regional scale because the 12th degree spherical harmonic solution is 
not adequate to establish 10® regional values. However, their spatial spectra do 
compare favorably (out to degree 12) since 10® regional values yield a reliable spec- 
tra to degree 12. Even though the shape factor and the deconvolution solutions com- 
pare most favorably on a 10® regional scale, the parameter estimation technique may 
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give more reliable 10® regional values than either of the two smoothed solutions 
However, because the true radiation field is unknown, one cannot demonstrate the 
accuracy of any of the three analysis techniques on the basis of these results. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 9, 1981 



APPENDIX A 


SPHERICAL HARMONIC COEFFICIENTS FROM REGIONAL VALUES 

Consider the radiant exitance field over the globe in terms of regional values 
which are the average radiant exitances over the regions. A set of regional values 
can be transformed to spherical harmonic coefficients by the principle of orthogonal 
projection, that is , 


^ ^ /sphere “ 


(A1 ) 


where C^ are the spherical harmonic coefficients and X(0,‘l>) is known in terms of 
regional values. These coefficients, however, represent the averaged or smoothed 
field. To compute the original, or unsmoothed, coefficients from the regional 
values, C™ must be increased, or enhanced, to represent the original field. The 
smoothing parameter defined by Pellinen (1967) is used for this purpose. Fol- 

lowing his approach, define the original spherical harmonic representation as 


X(0,^) 


GO 


E 

n=0 


N^^c”' cos rti^ + sin m$) 
1^0 ^ ^ ^ 


P^'Ccos 0) 
n 


(A2) 


and the spherical harmonic representation of the constant regional values which 
requires an infinite number of terms can be defined as 


X(0,^) 



cos m® + sin 
n 



P^(cos 0) 
n 


(A3) 


where X(0,^) = Xj^ is constant for < 0 < 0j^2 and ^ Rapp (1977) 

has shown that it is theoretically e3q>edient to use a circular smoothing operator in 
place of the rectangular regional operator (eq. (30)). The angular radius of the 
circular region Y* is chosen so that the areas of the circular region and the rec- 
tangular region are the same. This relationship is given 


. fHL\ _ / e sin 
\2 ) (,471 ) 

where 0 is the angular dimension of the rectangular region. A 10® x 10® rectangu- 
lar region corresponds to a circular region with Y* = 5.63®. 

The circular smoothing operator is similar to the rectangular smoothing operator 
(eq. (30)) and is given by 
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APPENDIX A 


/ 2n Y* 

f X(0,$) sin Y dY d^ 

fS 


^ 3=0 •'y=0 

k r2Tl .Y* 


r [ sin Y dY dP 
3=0 -^Y=o 


(A4) 


where Y and P are the polar coordinates describing the circular region centered 

0 +0 $ +$ 
k1 k2 , - k1 k2 


at 0, = 
k 

relationship : 


and ^ ^ Pellinen presented the following 


It 2*Jt 


I sin irt^ I 

m * 

p(Y) P (cos 0) 

'Y=0 P=0 ^ ) cos I 


/ L 


sin Y <iP dY 


211 p”^(cos 0. ) 
n k 


^sin m^_ ^it 

■ )/ 
^cos Y=0 


p(Y) P°(cos Y) sin Y dY 
n 


(A5) 


where the weighting function P(Y) is taken as 


P(Y) = 


(0® < Y < Y*) 

(Otherwise) 


Substituting equation (A2) into equation (A4) yields 


211 Y* n 

I ^ cos sin n^) p”^(cos 0) sin Y dY dP 

*^ft=n*^v=n »^=rn m=n ^ \ ^ n / n 


^ P=0‘"Y=0 n=0 m=0 

k 


211(1 - cos Y*) 


Simplifying with equation (A5) gives 


k‘T, r p •'"fc" “s 
ifei ^ " ‘ 


s”^ sin ^ P^^fcos 0, \ 
n k/ n\ k/ 


(A6) 


where 
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1 - cos Y* 


rY* 0 

/ P (cos Y^ sin Y ^Y 


"y=o 


Using 


P (X) 



2n + 1L^ 


P _,,{x) 
n+1 


d 

dx 



we can express the smoothing parameter as 


P 


n 


1 1 

1 - cos Y* 2n + 1 



(cos Y*) 


P (cos 
n+1 



(n — 1/ •••) 


and p = 1* 

0 

Now in equation (A6) represents the average over a circular region centered 

at 0j^ and If this expression is used to approximate the average over a rec- 

tangular region, then we substitute Xj^ for X(0,^) in equation (A3) and hy 

comparing coefficients obtain C™ = Pn^n ^n ^ Pn^n * Thus, the original spher- 
ical harmonic representation denoted by and is smoothed by the parameter 

P^ when averaged over a region. To compute the original spherical harmonic coef- 
ficients and from the regional values, we first must calculate the spheri- 

n n 

cal harmonic coefficients and S® that represent the smooth field by equa- 
tion (Al ) and then divide them by p^^; that is, = P^^ and sJJ = p^^ 

is interesting to note that as a result of equation (a 5), p^ is a function of 
degree n but not a function of order m of the spherical harmonic. The smoothing 
parameter p^^ is presented in figure 3 for two values of y** 
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APPENDIX B 


10® X 10® REGIONAL VALUES FROM 10® x 20® REGIONAL VALUES 

In this appendix, a formula for estimating 10® x io® regional radiant exitances 
from 10® X 20® regional radiant exitances is derived on a zonal basis so that the 
10® X 20® regional estimate of one zone affects only the 10® x lo® regional estimates 
of that zone. The criterion for dividing the estimate is that the s\im of the slopes 
squared around the zone is a minimum. This yields the solution with lowest spatial 
frequency. Consider a zone composed of 10® x 20® regions, where the radiant exit- 
ances are denoted ty M^, M 2 , ..., M^. Further let us define the 10® x 10 ® radiant 
estimates ty ^1 V ^12' ^22' ^2' where corresponds to the western 

half of the ith region and M ^2 corresponds to the eastern half. Thus, 

= (1/2)(M^-| + following diagram is helpful: 


M. 




^ r 




“11 

“12 

*^21 

“22 


|-^ ->■ 


M 


n 



1 

.'jij 


M. 


The sum of the slopes squared S is 



+ 


- \2 

\ ^ . 


2 


Substituting the relation M ^^2 ^ ” ^i1 U^ves 


S 



2M. 


M 



+ 



2M 
n 

M 



2 


or in matrix form 




APPENDIX B 


where 


= [M. M_ ... M ] 

1 ^ n 

^ ^21 • • • 


0 

0 


P = 


0 

1 


1 

0 


0 

0 


0 

1 


0 

0 


0 

0 


1 

0 


Since we desire to minimize S with respect to X, 


dS 

dX 


12X 


12m + 2PX + 2p\ - 4p"^M = 0 


or 


(61 + P + p'^)X = (61 + 2p'^)M 


(B1 ) 


The solution of equation (B1 ) for X defines (i = 1, 2, . . . , n) , and 

M .^2 “ These are the 10® x 10® regional radiant exitances that minimize 

the sum of the slopes squared and give the smoothest curve through the data. Solving 
equation (Bl) for n = 4 gives 


”ii 


= “ (M. . + 6M. - . ) 

6 x -1 1 1+1 


(i = 1, 2, 3, 4) (B2) 


where Mq = and . For n = 5, 


M 


11 


- 29 <-«i-2 * ="l-1 + - “i.1 




(i - 1, 2, 


5) 


where = ^-1 ^+2 ~ ^2* this study we have approximated the solution 

to equation (Bl) by equation (B2) for all n. 
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APPENDIX C 


SPECTRUM OF A RANDOM FIELD 


The interpretation of a nearly constant spectrum as random noise is illustrated 
by examining the spectrum of a random field. Define the random field by 5® x 5® 
regional values Xj^z where X^^ is a normal random variable with mean zero and vari- 
ance ^ • Also define E[X^^Xj] = 6^0^, which is the regional equivalent of white 
noise. From appendix A the spherical harmonic coefficients for this random field are 


and 


“ sin” K f 1®'®) 

" a Region k ■=” 


K)‘ ■ (=%)■ i i V, I 

and the expected value of ^^nj 


Region k •^Region j 


L 


y"* (0,$) (0',$') dA dA' 

cn cn 


r ,ii 

/ ^ \2 K 

^5 

■ (■•'.) £ 


i 


.m 


-i2 


Region k 


Y" (0,^>) dA 
cn 


(Cl) 


We derive a similar expression for E[(S^)^] and define 


r ^ / 

/ n.\2 



1 



+ E 

(=n) 
— — 

J 


(C2) 


—2 

which is represented in figure 6 for a = 10 W-m . Although the spectrum for the 
random field gradually increases for n > 15, it supports the interpretation of a 
constant spectr\am as noise. Note that if the regions are sufficiently small so that 
Y^n assxjmed constant over the region, then 


f y“ {Q,^) dA 

/__ _ cn 

‘^Region k 


' ' “^Region k ' ' 


dA 


and equation (Cl) becomes 
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E 


/ ^ \2 K 

PJ = W S i 


Region k 



For small regions the smoothing parameter is approximately unity for all n. 

Also, if we assume all areas to be equal, that is, = A = 4it/K, and make use of 

the orthogonality condition (eq. (23)), then 


K) 




X 


Sphere 



dA 


K 


with a similar e3q>ression for E[(sJJ)^], the expected value of the degree variance 
for small regions is approximately given by 


Thus, the spectrum of a random field defined on a regional grid system is propor- 
tional to }j2n + 1 , and for a 5® x 5® regional grid system the spectra given by equa- 
tions (C2) and (C3) are indistinguishable. 
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